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^ ! Abstract 
m 

^' \ We present a superfluid hydrodynamic model for the increase in moment of 

^\j ' inertia, A/, of molecules rotating in liquid ^He. The static inhomogeneous 

<y ■ He density around each molecule (calculated using the Orsay-Paris liquid 



o 



4 



He density functional) is assumed to adiabatically follow the rotation of 



On \ the molecule. We find that the AI values created by the viscousless and 

^i^ \ irrotational flow are in good agreement with the observed increases for several 

td ; molecules [OCS, (HCN)2, HCCCN, and HCCCH3]. For HCN and HCCH, our 

^ ■ model substantially overestimates A/. This is likely to result from a (partial) 

1-^ ■ breakdown of the adiabatic following approximation. 
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The spectroscopy of atoms and molecules dissolved in He nanodroplets provides both a 
new way to study microscopic dynamics of this unique quantum fluid [Q , and a very cold 
matrix (0.4 K 0) to create and study novel species [^^. Recent experiments have demon- 
strated that even heavy and anisotropic molecules display rotationally resolved vibrational 
spectra with a structure reflecting the gas phase symmetry of the molecule. However, the 
rotational constants required to reproduce the spectra are often substantially reduced from 
those of the isolated molecule. For example, the u^ vibrational band of SFg dissolved in He 
nanodroplets (first observed by Goyal et al. [^] and later rotationally resolved and analyzed 
by Hartmann et al. [0,[7[) indicates that the effective moment of inertia, /eg, in liquid ^He is 
2.8 times that of the isolated molecule. The same qualitative behavior has been found for a 
wide range of other molecules ||l]J^J^. In an elegant recent experiment, it has been demon- 
strated that the rotational structure of OCS broadens and collapses in pure ^He droplets, 
and is recovered when k. 60 "^He atoms are co-dissolved in the '^He [|l|. The association 
of the weakly damped, unhindered rotation with the Bose symmetry of ^He suggests that 
this phenomenon is a manifestation of superfluidity, and has been called the microscopic 
Andronikashvili experiment []I|. 

A theory able to reproduce the observed increase, AJ, in molecular moments of inertia 
would be of interest for at least two reasons. First, the enhanced inertia provides a window 
into the dynamics of the liquid. Second, the ability to predict the rotational constants would 
further improve the utility of He nanodroplet isolation spectroscopy for the characterization 
of novel chemical species. 

The first model proposed to explain the observed A/ assumed that a certain number of 
He atoms, trapped in the interaction potential of the solute, rotate rigidly with the latter . 
In the case of SFg, 8 He atoms trapped in the octahedral global potential minima would 
create a rigidly rotating 'supermolecule' that would have approximately the observed Jeg. 
In the case of OCS, putting a six He atom 'donut' in the potential well around the molecule 



also reproduces the observed Jeff |Tn]. Recent Diffusion Monte Carlo (DMC) calculations |TI 



have predicted that the effective rotational constant of SFg-HcAr monotonically decreases 
from that of the isolated molecule to the large cluster limit, reached at A^ = 8, and remains 
essentially constant for A^ = 8-20. The supermolecule model has been recently extended to 
consider the rigid rotation of a 'normal fluid fraction' of the He density which is claimed 
to be significant only in the first solvation layer p|,p!0|, based on Path Integral Monte Carlo 



calculations of Kwon et al. |T2[ which show a molecule-induced reduction of the superfluid 
fraction. These calculations have been recently used to propose a definition of a spatially 
dependent normal fluid fraction which reproduces the observed I^a of solvated SFg [l^ . 



The limitations of the supermolecule model are made clear by the A/ observed for HCN 
in He droplets p, a] which is only k, 5% of the A/ observed upon formation of a gas 
phase He-HCN van der Waals complex [Q. Furthermore, it has been previously recognized 
that in principle there is also a superfluid hydrodynamic contribution, /h, to A/ [|I],|15|,|16 



Previous estimates, based upon a classic treatment of the rotation of an ellipsoid in a fluid 
of uniform density, found that Jh is only a small fraction of the observed A/, at least for 



heavy rotors such as OCS |jTO|JT^. In this report, we show that if the spatial variation of the 
He solvation density around the solute molecule is taken into account, the calculated /h is 
instead rather large and agrees well with experimental data. We compare our calculations 
with the experimental results available in the literature (OCS [jT|, HCN R a]) and with 



results recently obtained in our laboratory for HCCCH3 and HCCCN, and in the laboratory 
of R. E. Miller for (HCN)2 g b] and HCCH |, c]. 

We first calculate the ground state He density, p, around a static solute molecule. The 
molecule is then considered to undergo classical rotation, slowly enough that the helium 
ground state density adiabatically follows the molecular rotation. The kinetic energy asso- 
ciated with the He flow (assumed viscousless and irrotational) is used to calculate I^. 

The main input of our hydrodynamic model is the ground state density of He around the 
solute. DMC calculations can provide this density with a minimum of assumptions beyond 
the interaction potentials [|r^, but are computationally expensive. The Density Functional 
method, which is a good compromise between accuracy and computational cost, consists in 
numerically minimizing the total energy of the many-body system in the form of a semi- 
empirical functional of the He density: E = J drT-C[p{r)]. The energy density Ti contains an 
effective non-local interaction with a few parameters fixed to reproduce known properties 
of bulk liquid He. The functional used here is the one termed Orsay-Paris |TB|, which was 



shown to accurately reproduce the static properties of pure and doped He clusters [jT9 
The need to treat axially symmetric molecules implies moving from one to two-dimensional 
equations. The new routines have been extensively tested against previously calculated 
spherically symmetric systems. The minimization of energy is carried out by mapping the 
density distribution on a grid of points and propagating it in imaginary time, starting from 
a trial distribution. 

The density functional also contains the interaction between the He and the impurity 
molecule. The interaction, assumed pairwise, is treated as a static external potential, since 
the molecules considered here are expected to have negligible zero point motion. Existing 



potentials for He-HCN [gg], He-HCCH |2T|, and He-OCS |2[ have been used without modi- 
fications. The He-(HCN)2 potential was generated as the superposition of the potential due 
to two HCN molecules whose centers of mass are separated by 4.44 A (the equilibrium dis- 



tance for the gas phase dimer PSf). The repulsive part of the He-HCCCN potential has been 



taken from p^; the attractive part from the He-HCN and He-HCCH potentials, using the 



concepts of distributed interaction and transferability |25]. The He-HCCH |21] and He-CH4 



26| interactions were used to generate the potential between He and HCCCH3, treating the 



latter molecule as cylindrically symmetric. Full detail on all potentials used are available 
from the authors, and will be published separately ^7 . 



Once the helium density profiles are calculated, the molecules are assumed to rotate 
perpendicularly to their symmetry axis with angular velocity u. We assume that the He 
density adiabatically follows this rotation, which allows us to calculate the laboratory-frame 
time-dependent density at each point in the liquid. This assumption is only valid if at each 
point the velocity of the fluid, v{r), is less than a critical velocity, Vc- If Vc is taken to be the 
velocity of sound, this is true for all our molecules, at the temperature of the droplet: 0.4 
K. A further justification to our assumption is also the fact that no critical value of angular 
momentum is experimentally observed for a wide class of molecules (i.e. for a wide range of 
fiuid velocities). 

The second assumption that we make is that the He behaves entirely as a superfiuid 
undergoing irrotational fiow. The assumption that the motion is irrotational implies that 
v(r) can be written as the gradient of a scalar potential: v = — V0 (the dependence of 
p, V, on r will be implicit from now on), where is known as the velocity potential. These 



assumptions lead to the following hydrodynamic equation for the velocity potential [p^] : 

V-(pV0) = |^ = -(Vp).(a;xr). (1) 

The first equality is just the continuity equation, while the second reflects the statement 
that the density is time-independent in the rotating frame. We select our axis system with 
z along the symmetry axis of the molecule, and assume that rotation takes place round 
the X axis with angular velocity a; = u; x. In order to better exploit the symmetry of the 
problem, we have used elliptical coordinates ^, Q, ip, where x = f ^/^'^^^sm{9) cos(v?), y = 
f \/^'^"-^sin(6') sin(y9), and z = f ^ cos{9). The surfaces of constant ^ are ellipses of rotation 
with foci at z = ±/. Two such surfaces limit the region where Eq. (|I|) is solved. The inner 
boundary excludes the volume occupied by the impurity, and is chosen as the largest ellipse 
contained in the region where p < O.OOSpo (po = 0.0218A~^ is the bulk liquid density). Von 
Neumann boundary conditions n ■ V0 = —n ■ {u x r) insure that the normal component 
of velocity matches the normal component of motion of the boundary [^. For the outer 



boundary, any ellipse large enough that the motion of the outside fluid is negligible can be 
chosen (with Dirichlet boundary conditions = 0). These boundary conditions result in a 
unique solution to the hydrodynamic equations. Other solutions exist if we do not require 
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the fluid to be irrotational, but it is known that these are higher in energy ||29|, and will 
include any solutions that have some portion (a "normal component" or a He "snowball") 
of the He density that rigidly rotates with the molecule. 

Given the solution, 0, to the hydrodynamic equation, we can calculate the kinetic energy, 
i^h; in the motion of the fluid by the following: 

(2) 
(3) 

Eq. @ follows directly from the definition of kinetic energy; Eq. (|^) is derived from Eq. (^ 
using standard vector identities and assuming that is a solution of Eq. (|I]). dS is defined as 
positive when pointing out of the region of the fluid, /h is the hydrodynamic contribution to 
the moment of inertia for rotation about the x axis, and mne is the atomic mass of helium. 
Both (dp/dt) and (p, are proportional to u, thus the above definition of I^ is independent of 
uj. The total kinetic energy of rotation will include the contribution from the molecule, K^ = 
^ImUj"^, where /m is the moment of inertia of the free molecule. We can also calculate the 
net angular momentum created by the motion of the He fluid: Jh = mne / pr x (— V0) dV. 
By use of standard vector identities and Eq. (^, this definition can be shown to lead to 
Jh = Ih<^- The total angular momentum is the sum of that of the rotating molecule and 
the total moment of inetia the sum of the moment of inertia of the molecule and that due 
to hydrodynamic motion of the superfluid. The local shape of the velocity field v(r) can be 
rather complex due to the presence of strong inhomogeneities in the density distribution. 

We calculate Jh by solving the hydrodynamic equation, Eq. (0) for 0, assuming unit 
angular velocity rotation around the x axis. It is computationally convenient to solve a 
slightly transformed version of Eq. (|I]), where the smoother function Inp(r) appears instead 
of p(r): 



VV+(Vlnp) ■ (V(/) + x X r) = 0. (4) 

Eq. (^) is solved, subject to the boundary conditions, by converting it to a set of finite 
difference equations on a grid of points in our elliptical coordinate system and using the 
Gauss-Seidel relaxatioia method ||30|. Both Eq. (|) and Eq. (H) are then evaluated by simple 
numerical quadrature, and are found to give the same value of Jh within a few percent. We 
also carefully tested the convergence of Jh with grid size. 

As an example of density distribution and velocity field. Fig. |1| shows our results for the 
OCS molecule in a cluster of 300 He atoms. On the left we give the contour plot of the He 
density near the molecule. One clearly sees the complex structure which results from the 
tendency to have He atoms near local minima of the impurity-He potential. The highest 
peak, at {y,z) = (—3.6, —1.2), corresponds to a ring of atoms perpendicular to the axis of 
the molecule. The integral of the density within this structure gives 6.5 atoms, and indeed 
7-8 is the number of He atoms one expects to fit into such a ring by close-packing. On 
the right side of the same figure we plot the current density, pv. We find that most of 
the kinetic energy density, |pv^, comes from the first solvation layer, the outer part of the 
cluster giving a negligible effect. 

In Table || our results are compared with existing experimental values for several 
molecules in He nanodroplets. There is an overall good agreement between the predicted and 
observed enhancements of the effective moment of inertia. From a quantitative viewpoint, 
one notices that the predicted moments of inertia tend to overestimate the experimental val- 
ues. In the case of the lightest rotors (HON and HCCH) the large discrepancy suggests the 
breakdown of the assumption of adiabatic following as recently predicted |ll|] . In that paper 



the importance of He exchange is pointed out; it is also shown that the interplay of the rota- 
tional constant with the potential anisotropy determines the extent to which the anisotropic 
He solvation density can adiabatically follow the rotation of the molecule. When the rota- 
tional constant of SFg is arbitrarily increased in the calculation by a factor of 10, the He 
density in the molecule-fixed frame becomes much more isotropic and the solvation-induced 



AJ decreases by a factor of 20 |Tl|]. We have recently obtained experimental evidence that 
AJ is larger for DON than for HON, which we believe to be direct experimental evidence 
for this effect 0. It is interesting to note that for these light (i.e. fast spinning) rotors the 
maximum of v(r) approaches the bulk ^He sound velocity. 

The overestimate of the moments of inertia for the other molecules likely reflects the 
uncertainties in the calculated p(r). We should remark here that while, by construction, 
the Orsay-Paris functional prevents p (the density averaged over an atomic volume) from 
becoming much larger than po, the functional was not constructed to deal with density 
gradients as high as those found in the first solvation layer. We observed that small changes 
in the form of the He density within the deep potential well of those molecules produce 
significant variations of the predicted moments of inertia, limiting the accuracy of the final 
results to 20% - 30%. This uncertainty does not affect the main result emerging from Table | 
that the hydrodynamic contribution to the moment of inertia of these systems, instead of 
being negligible, is rather large and can explain the observed rotational constants. 

One could object that the density values found at the minima of the He-molecule inter- 
action potential (e.g. ~ 11 po for OCS) are too high to be treated as those of a liquid, 
and should be interpreted as localized He atoms rigidly rotating with the molecule; it has 
been proposed that the He density distribution around the OCS-Heg supermolecule is only 



weakly anisotropic and thus can rotate without generating a significant hydrodynamic con- 



tribution []10l. We have calculated the above density distribution, and found that it is still 
strongly anisotropic, leading to a hydrodynamic moment of inertia of over 400 u ■ A . When 
combined with the moment of inertia of the OCS-Heg supermolecule, this gives a total ef- 

o 2 

fective moment of inertia of over 650 u • A , dramatically larger than the experimental value 
(230 u- A'). 

In summary, the spatial dependence of the He density, which is caused by the molecule- 
He interaction, results in a hydrodynamic contribution to the moment of inertia more than 
an order of magnitude larger (in the case of the heavier rotors) than that predicted for the 
rotation of a reasonably sized ellipsoid in He of uniform bulk liquid density. Furthermore, 
the present calculations suggest that the effective moments of inertia of molecules in He nan- 
odroplets (and likely also bulk He) can be quantitatively predicted by assuming irrotational 
fiow of a spatially inhomogeneous superfiuid. 
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REFERENCES 

[1] S. Grebenev, J. P. Toennies, and A. F. Vilesov, Science 279, 2083 (1998). 

[2] M. Hartmann, R. E. Miller, J. P. Toennies, and A. Vilesov, Pliys. Rev. Lett. 75, 1566 
(1995). 

[3] J. Higgins et al, Science 273, 629 (1996). 

[4] K. K. Lehmann and G. Scoles, Science 279, 2065 (1998). 

[5] K. Nauta and R. E. Miller, Science 283, 1895 (1999). 

[6] S. Goyal, D. L. Schutt, and G. Scoles, Phys. Rev. Lett. 69, 933 (1992). 

[7] R. Frochtenicht, J. P. Toennies, and A. F. Vilesov, Chem. Phys. Lett. 229, 1 (1994). 

[8] K. Nauta and R. E. Miller, (a) Phys. Rev. Lett. 82, 4480 (1999); (b) J. Chem. Phys. 
Ill, 3426 (1999); (c) private communication. 

[9] A. Conjusteau et al, in preparation. 
[10] S. Grebenev et al, preprint. 

[11] E. Lee, D. Farrelly, and K. B. Whaley, Phys. Rev. Lett. 83, 3812 (1999). 
[12] Y. Kwon, D. M. Ceperley, and K. B. Whaley, J. Chem. Phys. 104, 2341 (1996). 
[13] Y. Kwon and K. B. Whaley, Phys. Rev. Lett. 83, 4108 (1999). 
[14] S. Drucker, F. M. Tao, and W. Klemperer, J. Phys. Chem. 99, 2646 (1995). 
[15] K. B. Whaley, Advances in Molecular Vibrations and Collision Dynamics, eds. J. M. 

Bowman and Z. Bade, 397 (JAl Press, 1998). 
[16] K. K. Lehmann, Mol. Phys. 97, 645 (1999). 

[17] R. N. Barnett and K. B. Whaley, J. Chem. Phys. 99, 9730 (1990). 
[18] J. Dupont-Roc, M. Himbert, N. Pavloff and J. Treiner, J. Low Temp. Phys. 81, 31 

(1990). 
[19] F. Dalfovo, Z. Phys. D 29, 61 (1994). 

[20] K. M. Atkins and J. M. Hutson, J. Chem. Phys. 105, 440 (1996). 
[21] R. Moszynsky, P. E. S. Wormer, and A. van der Avoird, J. Chem. Phys. 102, 8385 

(1995). 
[22] K. Higgins and W. Klemperer, J. Chem. Phys. 110, 1383 (1999). 
[23] K. W. Jucks and R. E. Miller, J. Chem. Phys. 88, 6059 (1988). 
[24] E. M. Cabaleiro-Lago and M. A. Rios, J. Chem. Phys. 109, 8398 (1998). 
[25] R. J. Bemish, R. E. Miller, X. Yang, and G. Scoles, J. Chem. Phys. 105, 10171 (1996). 
[26] U. Buck et al, Mol. Phys. 55, 1255 (1985). 
[27] C. Callegari et. al, in preparation. 

[28] L. M. Milne-Thomson, Theoretical Hydrodynamics, fifth ed. (Dover Publications, Inc., 
Mineola, NY, 1996). 



[29] Kelvin's minimum energy theorem. The proof in: [28] is easily extended to the non- 
uniform density case. 
[30] Numerical Recipes, Chapter 19, Cambridge University Press. 



TABLES 





/m 


Ics (exp.) 


A/ 


Ih (calc.) 


ref. 


HCN 


11.39 


14.04 


2.65 


5.47 


I| a] 


DCN 


14.0 


16.9 


2.9 


5.6 


|27| 


HCCH 


14.26 


16.08 


1.82 


9.4 


II c] 


HCCCH3 


59.14 


224.0 


164.8 


190 


11 


OCS 


83.10 


230.0 


146.9 


197 


II 


HCCCN 


111.1 


330.7 


219.6 


226 


11 


(HCN)2 


289.5 


872.5 


583 


619 


II b] 
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TABLE L Moments of inertia for the molecules studied in this work. Units are u • A . The 
quantities /m and leflf are the observed moments of inertia when the molecule is free and dis- 
solved in the cluster, respectively. Their difference, A/, in the 4*^ column is compared with the 
hydrodynamic moment of inertia, /h; of the present calculation. 



FIGURES 

FIG. 1. He density, p, (left) and He current density, pv, (right) distributions for a cluster of 
300 He atoms, with OCS in its center rotating counterclockwise. For the sake of showing details, 
the highest He density peaks have been clipped (white areas), and the dynamic range of the current 
has been compressed. 



